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We investigate spin properties of a Fe^”*" dopant, known for having single nondegenerate ground 
state in bulk host semiconductor. Due to zero magnetic moment such a ground state is of little use 
for spintronics and solotronics. We show that this well-established picture of Fe^”'" spin conhguration 
can be contradicted by subjecting the Fe^"*" ion to sufficiently high strain, e.g., resulting from lattice 
mismatched epitaxial heterostructures. Our analysis reveals that high strain induces qualitative 
change in the ion energy spectrum and results in doubly degenerate ground state with spin projection 
Sz = ±2. An experimental proof of this concept is demonstrated using a new system: an epitaxial 
quantum dot containing individual Fe^^ ion. Magnetic character of the Fe^"*" ground state in a 
CdSe/ZnSe dot is revealed in photoluminescence experiments by exploiting a coupling between a 
confined exciton and the single iron impurity. 

PACS numbers: 78.67.He, 78.55.Et, 75.75.-c 75.30.Hx 


Spin configurations of transition metal ions in various 
host semiconductors have been well established already 
a few decades ago [H-Q . It has been found that ions such 
as Gr^+((i'^), Mn^+(d®), Co^~^{d7) exhibit nonzero spin 
in their ground states, which makes them useful in spin¬ 
tronics [3 ,[i 3- However, the ground state of the Fe^'''(d®) 
ion in zinc-blende or wurtzite II-VI semiconductors like 
ZnS, ZnSe, GdTe or CdSe has been found to be nonde¬ 
generate 11 and thus termed nonmagnetic (^ . To 
induce Fe^“'" magnetic moment, high magnetic field has 
to be applied, as for Van Vleck paramagnets 21, 

The physics of the transition metal ions has been re¬ 
cently brought back into the spotlight due to possibil - 
ity to access to the properties of single dopants 
Among other achievements, optical orientation 
readout 24, M [s^l and coherent precession [S^ of 
a single magnetic ion spin were demonstrated. Current 
development of the field benefits greatly from the fun¬ 
daments of the early findings. However, the physics of 
the transition metal ions in semiconductor nanostruc¬ 
tures goes far beyond the limits established in the earlier 
works on bulk materials. 



In this Letter we demonstrate that, contrary to the 
well-established knowledge on a Fe^'*' ion in the semi¬ 
conductor matrix, it is possible to qualitatively change 
its low-field behavior from nonmagnetic to magnetic, in 
particular by placing such an ion in a highly strained 
nanostructure. In order to elucidate this fact, we analyze 
the Fe^"*" energy spectra for the cases of weak and strong 
strain, showing a clear hierarchy of the energy scales, rel¬ 
evant both to zinc-blende and wurtzite structures. The 
magnetic behavior of the Fe^“'" ion is experimentally ev¬ 
idenced by analyzing the magnetic field dependence of 
the photoluminescence (PL) spectrum of an individual 
GdSe/ZnSe quantum dot (QD) containing a single Fe^+ 
impurity. The nonzero spin ground state of the Fe^“'" ion 
opens the possibility of using it as a two-level system in 


quantum information technology [3- 


The dominant effect defining energy spectrum of a 
transition metal ion in the bulk semiconductor is the crys¬ 
tal field [1, [l3|- Fe^'*' has configuration d®, which means 
that the d-shell electrons have combined orbital angular 
momentum of L = 2 and spin of S' = 2. The crystal field 
affects only orbital part of the wave function and splits 
five orbital states of the ion into two subspaces: twofold 
degenerate and threefold degenerate ^T 2 , with be¬ 
ing lower-energy in Td symmetry (Fig. [Ija)). For Fe^“'" 
in CdSe or ZnSe this splitting is about 10\Dq\ fv 0.3 eV 
SElsi. Thus, the ^T 2 level is not populated even at 
room temperature and the properties of the Fe^+ ion are 
defined only by the states in the subspace. These 
states are not affected by a static Jahn-Teller distor- 
tion, as it was shown for many Fe-doped semiconductors 
[si IiiI413| . Gonsequently, the second effect in order of 


strength is the spin-orbit interaction XLS. It results in 
splitting of ® A levels into 5 equidistant groups, as shown 
in Fig. HJa). The value of the splitting is given by the 
effective strength of XLS interaction and the crystal field 
splitting: ATls = GX"^/10\Dq\ Ri 2 meV [l3 O- The 
presence of a dynamical Jahn-Teller effect or application 
of additional stress in experimentally accessible range re¬ 
sults in only small shifts of those energy levels and can 
be treated perturbatively 0,0, la every case, 

the lowest-energy group consists of a single nondegener¬ 
ate state, which determines the nonmagnetic character 
of the Fe^’^' ion ground state. 


We find that strong structural strain of a QD changes 
hierarchy of the Fe^'*' energy scales. The dominant effect 
is still the crystal field, but the second effect becomes the 
biaxial strain. It lifts orbital degeneracy of the sub¬ 
space, splitting it into states of symmetries correspond¬ 
ing to single-electron and orbitals (Fig. HKb)). 

The ordering of those states is determined by the sign 
of the strain, which, given the GdSe/ZnSe lattice mis- 
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FIG. 1: (Color online) Energy spectrum of a Fe^^ ion in (a) 
a bulk zinc-blende semiconductor, (b) a nanostructure with 
a strong in-plane compressive strain. Numbers in parenthe¬ 
ses denote the degeneracy of the energy levels. Labels of the 
orbital states split by the QD strain refer to single-electron 
orbitals of corresponding symmetry (for more details, see Sup¬ 
plemental Material), (c) Visualization of two orbitals forming 
the subspace. Arrows schematically mark the shift of the 
neighboring anions due to the strain of the QD. 


match, has compressive character. Qualitatively, such 
strain pulls the tetrahedral lattice bonds away from xy 
plane and thus lowers the energy of the orbital 

while increasing the energy of the (ij ,2 one (as schemati¬ 
cally depicted in Fig. HKc)). More strict analysis leading 
to the same level ordering is presented in Supplemental 
Material (SM). Finally, the spin part of the wave function 
is determined by the XLS interaction. This interaction 
contributes to the energy of spin states of dx 2 _y 2 orbital 
in the second order. It favors high spin states according 
to the effective spin Hamiltonian DS^ with D < 0 (for 
details, see SM). As a result, the ground state is dou¬ 
bly degenerate (within the discussed second order of the 
XLS interaction) with the spin part of Sz = ±2. Such 
two states are easily split by external magnetic field, in 
a clear contrast to previously described case of the Fe^+ 
embedded in bulk semiconductor. 

The experiment evidencing actual behavior of the Fe^“'" 
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FIG. 2: (Golor online) (a) A PL spectrum of a typical CdSe 
QD showing neutral exciton (X), charged exciton (X~), and 
biexciton (2X) lines. Neutral complexes exhibit anisotropic 
splitting of (5i = 370 yeV. (b) A PL spectrum of a QD with 
a single Fe^”*" ion. The PL lines are split mainly due to s,p-d 
exchange interaction between confined carriers and the d-shell 
electrons of the ion. For both spectra continuous background 
was subtracted. Inset: histogram of the s,p-d exchange split¬ 
ting of the X“ emission line. The cut-off at As,p-d 0.3 meV 
is due to our selection procedure related to the resolution 
of our experimental setup — only dots with larger zero-field 
splitting were tested in the magnetic field to verify the pres¬ 
ence of the Fe^^ ion. 


ion in a strained nanostructure is carried out on a number 
(> 30) of single QDs, each incorporating an individual 
iron ion. Such structures are presented here for the first 
time. Self-assembled zinc-blende CdSe QDs in ZnSe bar¬ 
rier are grown using molecular beam epitaxy. About 2 
monolayers of CdSe:Fe are deposited without any growth 
interruptions on 1.5 ym thick ZnSe buffer grown on GaAs 
(100) substrate. QDs are covered by a 50 nm thick ZnSe 
cap layer. Iron doping density is adjusted in order to op¬ 
timize the probability of finding a QD with exactly one 
Fe^'*' ion. Low temperature (^1.5 K) PL experiments on 
individual QDs are performed in a setup providing spa¬ 
tial resolution of 0.5 ^m without the need for mesas or 
masks. The PL is excited nonresonantly at 405 nm. The 
magnetic field up to 10 T is applied in Faraday configu¬ 
ration. 

As expected for random character of low density dop¬ 
ing, in the same sample we find QDs incorporating single 
Fe2+ ions and undoped QDs for reference purposes. PL 
spectra corresponding to both of these cases are shown 
in Fig. [2] Figure [2Ka) presents a typical spectrum of 
an undoped QD. The spectrum exhibits all standard fea¬ 
tures of epitaxial quantum dots [H, IdU - li^ l . The sharp 
emission lines originate from recombination of differ- 
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ent excitonic complexes, including neutral exciton (X), 
negatively charged exciton (X“), and biexciton (2X). 
The neutral exciton and biexciton lines are split due to 
anisotropic part of electron-hole exchange interaction. In 
the case of QD shown in Fig. HJa) this splitting yields 
(5i = 370 /reV, which is a representative value for our sam¬ 
ples. On the other hand, the charged exciton line does 
not exhibit any splitting, in accordance with the Kram- 
mers rule for systems with odd number of fermions. 

In comparison, introduction of a single Fe^'*" ion into a 
QD leads to distinctive changes in the PL spectrum, as 
shown in Fig. [2l[b). The emission lines still correspond 
to recombination of the same excitonic complexes, how¬ 
ever their structure is determined by the s,p-d exchange 
interaction with the resident ion. The main effect is a 
strong splitting of each of the observed emission lines. It 
is particularly striking for typically degenerate charged 
exciton, but also for the neutral exciton it is significantly 
stronger than typical value of . Such a physical picture 
is similar for large number of studied Fe-doped QDs, as 
proven by distribution of measured s,p-c? exchange split¬ 
tings presented in the inset to Fig. [2jb). The presence of 
such s,p-d splitting is a direct confirmation of the mag¬ 
netic character of the Fe^+ ion. It originates from the 
fact that the Fe^"*" spin may be aligned either parallel 
or anti-parallel to the exciton angular momentum, which 
would not be possible in the case of nonmagnetic ground 
state. 

In order to provide the final proof of the magnetic char¬ 
acter of the Fe^+ ion in a QD, we measure the evolution 
of the X PL spectrum in external magnetic field applied 
along the growth direction (quantization axis of the mag¬ 
netic ion and QD excitons). Typical results obtained in 
a~ polarization of detection are shown in Fig. |3l(a). We 
note that the observed pattern is quite similar to the 
one obtained for InAs/GaAs QDs containing single man¬ 
ganese ions [ 23 , Hi ) despite different microscopic origin. 

Magneto-photoluminescence results in Fig. [31(a) seem 
complex, however they originate from quite simple be¬ 
havior of the initial and final energy levels of the tran¬ 
sitions, as illustrated in Fig. |3 Kc). First effect of the 
magnetic field is the Zeeman splitting between Sz = 2 
and Sz = —2 states of the Fe^’*' ion. Unfortunately, the 
PL spectrum does not show this splitting directly, since 
in general exciton recombination does not influence the 
ion spin state and thus the energy of emitted photon does 
not depend on the ion Zeeman splitting. However, in the 
vicinity of the level anticrossings the Fe^+ spin states are 
mixed and this selection rule is relaxed. Indeed, data in 
Fig. ETa) features several weaker lines in the anticrossing 
range (i.e., 0 - 2 T). Before we discuss the origin of the 
anticrossings, let us analyze the behavior of these weak 
lines, in particular the cross-like feature. The two cross¬ 
ing lines correspond to transitions involving the change of 
the ion spin from Sz = ±2 to Sz = t 2. The splitting be¬ 
tween them depends almost linearly on the magnetic field 


FIG. 3: (Color online) Magnetic field dependence of the PL 
spectrum of a neutral exciton in a QD with a single Fe^”'" ion: 
(a) experimental data, and (b) simulation assuming strain- 
induced magnetism of the Fe^"*" ion, as described in the text. 
The spectra were measured (simulated) in <j~ circular po¬ 
larization. (c) Schematic dependence of the involved energy 
levels on the magnetic field together with indicated a~ polar¬ 
ized optical transitions observed in PL measurements. The 
upper pair of levels corresponds to IL-II) exciton coupled with 
the ion spin (where L and JJ. represent the spin projection of 
the electron and the heavy hole on the growth axis, respec¬ 
tively), while the bottom pair represents the energies of the 
ion states in the empty dot. The excitonic transitions preserv¬ 
ing (altering) the ion spin projection are marked with solid 
(dashed) arrows. 


with a slope of about 0.84 meV/T. More precise fitting 
including non-linearity due to proximity of the anticross¬ 
ings gives slightly larger value of 0.92 meV/T. Taking into 
account that 1X52 1 = 4 for both the initial and the final 
states, this slope corresponds to g-factor of 2.0, exactly 
as expected for the Fe^+ spin. 

Let us now focus on the nature of the observed anti¬ 
crossings. The first, relatively weak anticrossing occurs 
around 0 T. It is a signature that the Sz = ±2 states of 
the Fe^’*' ion are not perfectly degenerate, but are split 
by a small energy o, as shown in Fig. EIc). This split- 
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ting varies between different studied dots and its average 
value yields about 50 /xeV. Consequently, the presence 
of this splitting does not invalidate our conclusion about 
magnetism of the ground state of the Fe^+ ion, since a 
is negligible even compared with the X-Fe^'*' exchange. 
Such a splitting requires including XLS coupling in the 
fourth order of perturbation theory, according to which 
the zero-field eigenstates are -^{\Sz = 2) ± = — 2)). 

The second anticrossing around 2 T is closely related 
to the first one. It occurs when the effective magnetic 
field acting on the Fe^’*' spin in the presence of the a~ 
emitting exciton vanishes. Since exchange field of this 
exciton increases the energy of the state corresponding to 
Sz = —2 ion spin projection (Fig. [3Kc)), the anticrossing 
of the Fe^"*" ion is effectively shifted from 0 T to a higher 
field. 

Finally, there is also the third, stronger anticrossing 
around ±9 T. This anticrossing is observed for both neg¬ 
ative and positive magnetic field (or equivalently: in a'^ 
and a~ polarization), which clearly indicates that it is 
due to mixing of the exciton part of the total wave func¬ 
tion. Indeed, the states involved in the anticrossing cor¬ 
respond to cr~ and emitting excitons coupled with 
Sz = —2 spin projection of the Fe^+ ion (ItJJ-,—2) and 
|4,jl, —2)). The anticrossing occurs when the excitonic 
Zeeman effect reduces the ion-related exchange splitting 
of the involved states and the anisotropic electron-hole 
exchange interaction becomes dominant source of the 
splitting. It should be noted that this anticrossing does 
not mix different states of the Fe^"*" ion and thus in this 
range of magnetic held the optical recombination pre¬ 
serves the spin of the ion. 

In order to quantitatively verify our interpretation of 
the magneto-photoluminescence results, we perform a 
numerical simulation of the expected held-dependence of 
X PL spectrum. The simulation is based on the general 
spin Hamiltonian of an ion-exciton system 23|, l24l . 


H — Hio^ + Hx + Hg p-d+fiBB (^ion'S'z™ + geS^ + OhS^') , 


where Llion is the Hamiltonian of the Fe^’*' ion leading to 
energy spectrum as in Fig[T](b), Hx is the Hamiltonian of 
the exciton with electron-hole exchange interaction 4Ij . 
Hg^p-d describes the s,p-d exchange interaction between 
conhned carriers and the iron ion (dH - ldojl . and gion, 5e 
and gh are g-factors of the Fe^+ ion, electron and hole, 
respectively. As shown in Fig. [3Kb), such a model re¬ 
produces all features of the experimental measurement 
(for details of the simulation procedure, see SM). The 
model correctly captures even the observed thermalisa- 
tion of the ion spin at increasing magnetic field by taking 
into account the effective Fe^’*' spin temperature of 15 K. 
Such a good overall agreement provides a strong proof of 
correct identification of all relevant effects. 

All the presented results clearly show that the struc¬ 
tural strain of the QD induces magnetic character of the 


Fe2+ ion in its ground state. Our findings reveal that 
a CdSe/ZnSe QD containing such an ion may be use¬ 
ful for spin-based quantum information technology, as it 
combines many desired qualities. In particular, both the 
CdSe lattice and Fe^"*" ion can be free of any nuclear spin 
fluctuations. Moreover, the QD provides efficient opti¬ 
cal access to a single ion. Finally, as we show here, the 
ground state of the Fe^"*" ion has nonzero spin, which 
opens the possibility of using such ion as a two-level sys¬ 
tem. However, the importance of our results is not lim¬ 
ited to this particular system. It is a general example of 
the fact that even well-established textbook knowledge of 
energy spectrum of various dopants should be carefully 
re-evaluated in the world of semiconductor nanostruc¬ 
tures. 
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LABELING OF ION ORBITAL STATES 

In Table H] we list the symmetries of five multi-electron 
orbital states of Fe^'*' ion. Since six d-shell electrons have 
total orbital angular momentum L = 2, their orbital wave 
functions transform like linear combinations of spherical 
harmonics , where = 0, ±1, ±2. They are explic¬ 
itly given in the second column of Table U and denoted by 
symbols rj, C, 0, e after J. T. Vallin et al. Q. Since d® 
conhguration corresponds to half-filled shell with one ad¬ 
ditional electron, these multi-electron orbital states may 
be also discussed using single-electron d orbitals of the 
same symmetries, which are listed in the last column. 


Subspace 

Orbital wave function 

Corresponding 

single-electron 

orbital 

®T2 

|,7) = -(y,i-y-i)/V2 

dxz 

10 = i{Yi + Y^-^)/V2 

dyz 

10 = iiXi - 

dxy 


II 

d^2 

10 = {Yi + Y^-^)/^ 

dx2_y2 


TABLE I: Symmetry of Fe^'*' orbital states. 


ENERGY LEVELS OF FE2+ ION IN THE 
PRESENCE OF BIAXIAL STRAIN 


In order to estimate the magnitude of strain-induced 
splitting of orbitals we employ a simple point-charge 
model, in which the crystal field acting on the single 
Fe^"*" ion comes from four point charges corresponding to 
neighboring Se^“ anions (Fig. [Ha)). The electrostatic 
potential from these charges can be expanded in multi¬ 
pole series. In the case of pure tetrahedral Td symmetry, 
the terms relevant to d orbitals of the Fe^’*' ion are given 
by 


Vcir) = A 


d4 




( 1 ) 


where are spherical harmonics in a coordinate 

system (r, 6, tf) centered on the Fe^+ ion. The coefficient 
A denotes the strength of the crystal field and is given by 
A = d ^ where Cr is the relative permittivity 



FIG. 1: (a) A schematic visualisation of Fe^”*" ion snrrounded 
tetrahedrally by four Se^“ anions occupying alternating cor¬ 
ners of a cnbic elementary cell, (b) Illustration of strain- 
induced displacement of each Se^“ anion (not in scale). The 
angle 8 corresponds to the anion position in pnre Td symme¬ 
try, while A9 denotes the change arising due to the strain. 


of CdSe taken as 9.2 Q and d = a\/3j2 is the distance 
between Fe^+ and each anion in the crystal lattice (with 
a ~ 6.08 A being zinc-blende GdSe lattice constant Q). 

The structural strain of a CdSe/ZnSe QD shifts the 
anions away from their normal positions and thus lifts the 
tetrahedral symmetry. The main factor determining the 
character of the strain is a large lattice mismatch between 
CdSe and ZnSe, which is equal to about 7% [i,[3. As 
a result, the dot material is compressed in xy plane and 
stretched along the growth direction z. We assume that 
the relative reduction of a lattice constant is the same 
for X and y directions and equal to 5x- Gonsequently, 
the relative elongation of an elementary cell in z is given 
W = 2u5x, where v k, 0.35 is GdSe Poisson’s ratio 
[sj. Under these assumptions, the angular coordinate ip 
of each anion in xy plane remains unchanged, while the 
angle 9 between the z axis and direction of Fe^^-Se^" 
bond is slightly decreased, as schematically shown in Fig. 
[IJb). To the accuracy of linear terms in 5x, the change 
of 9 angle is given by A0 = \/2(l -I- 2v)5x/i, which for 
maximal 5x limited by the lattice mismatch corresponds 
to about 3°. 

Such strain-induced displacement of the anions lowers 
the symmetry to D 2 d and introduces new terms to the 
crystal field potential, which for d orbitals read 


Vs{r)=Ae~ 

0 _ 5^5 
^ 27 
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where Ag Ri 16 A0. The overall potential 

Vc + 14 ' splits fivefold degenerate orbitals of the Fe^"*" 
ion into four subspaces, which energies take a form 


E(dxz) = E{dyz) 

= —ADq + Ds 2Dt, 

(3) 

E{dxy) 

= —ADq — 2Ds — ADt, 

(4) 

E{dz2) 

= 6Dq + 2Ds — 3Dt, 

(5) 

E{dx'2 — y2 ) 

= 6Dq — 2Ds + 3Dt, 

(6) 


where 10|Zlq| (with Dq < 0) is a splitting between ®T 2 
and subspaces arising only due to the crystal field of 
pure Td symmetry, while Ds and Dt are the splittings 
induced by the strain. They are given by 


Dq = - 


Ds = 1^(1 
21^ 

160 


27 dTrepCrd^ 

2v)5: 




Dt = 


567 


(1 + 2v)5x 


iTTeoer(P 
^2 


Aireoerd^ 




(7) 

( 8 ) 

(9) 


where (r") denotes the expectation value of nth power 
of Fe^+ orbital radius r. According to Eqs. ([S]) 
and (0), the strain-induced energy splitting of orbitals 
dx^-y 2 and dz^ forming the subspace corresponds 
to As{^E) = ADs — 6Dt. Under an assumption that 
(r") R! (r)" it can be expressed using \Dq\ and reads 


Asi^E) = {1 +2,^)6, —\Dq\- 

(^ 

y 5 y 47reoerd|£>g| 



Since lO\Dq\ r 0.3 eV for CdSe [§-[lll, the splitting 
Asi^E) corresponds to 4 • 0.7 eV, which for maximal 6x 
of 0.07 yields about 50 meV. This value is significantly 
larger than the splitting oA^E orbitals arising from spin- 
orbit coupling with higher energy states, which was found 
to be of the order of a few meV in bulk semiconductors 
[lll - [l3j| . Consequently, our calculations reveal that the 
crystal field and strain are dominant effects influencing 
the energies of the Fe^"*" ion states in a CdSe/ZnSe QD 
and the spin-orbit interaction might be treated pertur- 
batively. Moreover, our results also confirm that d^'^-y^ 
is the lowest-energy orbital state of this ion embedded 
in a dot, since the computed Agi^E) is positive for the 
compressive strain (i.e., for 4 > 0). It should be noted 
that the actual value of this splitting may be deviated 
to a some degree from the one given by Eq. , since 
the exploited point-charge model neglects the covalency 
of bonds between magnetic ion and anions in the crystal 
lattice. Nevertheless, as proven by the experiment, the 
general conclusions of our analysis are not altered by this 
simplification, even though it is probably quite impor¬ 
tant, especially in the case of II-VI compounds studied 
in this work. 


EFFECTS OF SPIN-ORBIT COUPLING 

Since both the crystal field and strain does not af¬ 
fect the spin part of the Ee^'*' ion states, the lowest- 
energy orbital d^^-y^ is fivefold degenerate due to the 
ion spin S = 2. This degeneracy is partially lifted by the 
spin-orbit interaction XLS. Such interaction splits the 
spin states of d^'^-y^ orbital in the second order through 
mixing with higher energy orbitals forming the ^T 2 sub¬ 
space. This splitting can be expressed by the effective 
spin Hamiltonian: 


Hls = D 



^*5(5+1) 


( 11 ) 


where D is given by 


D = 


A2 

E{d 

XZ )-E{dx2.y2) 


4A2 

Ei^dxy^ E(^dx2—y2^ 


( 12 ) 


Taking into account Eqs. @-(15]) and bearing in mind 
that Ds and Dt are positive for the compressive strain, 
one obtains E{dxz) > E{dxy) > E{dx^-y' 2 ), which yields 
D < /[E{dxy) — E{dx 2 -y 2 )] < 0. Consequently, the 
states corresponding to Sz = ±2 spin projections have 
lowest energy. Their degeneracy is however lifted, when 
the spin-orbit coupling is considered in the higher order. 
In particular, in the fourth order the states |S '2 = ±2) 
are split by XLS interaction into linear combinations 
-^{\Sz = 2) ± \Sz = —2)), with the state corresponding 
to being lower energy. The corresponding energy 
splitting a is given by 


36A4 

[Eidxz) - Eidx2_y2)f Asi^E)' 


(13) 


The value of this splitting obtained for maximal strain 
(i.e., 4 = 0.07) and A r -10 meV [i [H is equal to 
about 55 /ieV, which is consistent with our experimental 
results. 


SIMULATION OF THE PL SPECTRUM OF A QD 
WITH A SINGLE FE2+ ION 

As shown in the manuscript, all effects observed in the 
magnetic field evolution of the PL spectrum of a QD with 
Fe^“'" can be reproduced by a numerical simulation. Our 
simulation is based on the standard procedure of find¬ 
ing the eigenstates of the exciton complex and the empty 
(i.e., without the exciton) QD and subsequent calculation 
of allowed transitions. For simplicity, we assume that the 
spatial part of exciton wavefunction is not substantially 
modified by the magnetic field and we can restrict our 
analysis to the spin degree of freedom. The possible ini¬ 
tial states were found by diagonalization of the Hamilto- 
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nian: 

Hi = -ffion + Hx + Hs^p-d + 

+ f^BB{gionS'r + geSl+gi,S^) (14) 

where S''™, S®, S*' are effective spin operators of the ion 
and carriers with eigenvalues of ±2, ±1/2, ±3/2, respec¬ 
tively. The Hamiltonian is diagonalized within the space 
of possible orientations of bright exciton and the Fe^’*' 
spin: |tJ/,±2), |tJ/,-2), lilt, ±2), lilt,-2). Dark exci- 
tons were not included in the simulation, since due to 
the axial symmetry of the Hamiltonian they do not mix 
with the bright excitons. 

The Hamiltonian of the final state includes only the 
contributions from the Fe^^ ion: 

Hf = FTion ± ffionMsSSr. (15) 

The Hamiltonians are parameterized by 5 free param¬ 
eters: di, gx = 5gh — ffe, 5ion, the zero-field splitting a of 
the Fe^’*' states, and the exciton-ion exchange splitting 
As^p-d.. The values of these parameters are found by fit¬ 
ting the field dependence of the PL spectrum. Additional 
parameter Eq is introduced later to plot the simulation 
in the same energy range as the experimental results. 

Since spatial part of all considered states is the same, 
the oscillator strength between an initial and a final state 
is determined only by the spin part of the wavefunction. 
In our case the oscillator strength is proportional to the 
projection of the initial state |f) onto the state corre¬ 
sponding to the final spin projection of the ion |/) (e.g., 

I±2)) and the given orientation of the exciton spin (e.g.. 
It'll) for simulation of a— spectrum). 

Finally, the optical spectrum is simulated as a series 
of peaks at energies corresponding to possible transitions 
between possible initial and final states. The intensity 
of a given peak is the product of calculated oscillator 
strength and the Boltzmann factor corresponding to ther- 
malisation of the ion spin at an effective temperature of 


15 K. For the sake of the presentation, the peaks are 
plotted as gaussians with FWHM of 0.1 meV. Strictly 
speaking, such procedure gives us the low-power absorp¬ 
tion spectrum, but it is also similar to the PL spectrum 
if only the exciton states have the same lifetime and that 
is shorter than the thermalisation of the excited state. 
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